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Abstract 

We investigate resonant spin-flavor (RSF) conversions of supernova neutrinos which are induced 
by the interaction of neutrino magnetic moment and supernova magnetic fields. From the formu¬ 
lation which includes all three-flavor neutrinos and anti-neutrinos, we give a new crossing diagram 
that includes not only ordinary MSW resonance but also magnetically-induced RSF effect. With the 
diagram, it is found that four conversions occur in supernova, two are induced by the RSF effect and 
two by the pure MSW. We also numerically calculate neutrino conversions in supernova matter, us¬ 
ing neutrino mixing parameters inferred from recent experimental results and a realistic supernova 
progenitor model. The results indicate that until 0.5 seconds after core bounce, the RSF-induced 
v e <-> v T transition occurs efficiently (adiabatic resonance), when p v > 10 -12 pb(Bq/5 x 10 9 G) -1 , 
where Bq is the strength of the magnetic field at the surface of iron core. We also evaluate the en¬ 
ergy spectrum as a function of g u Bo at the SuperKamiokande detector and the Sudbury Neutrino 
Observatory using the calculated conversion probabilities, and find that the spectral deformation 
might have possibility to provide useful information on neutrino magnetic moment as well as mag¬ 
netic field strength in supernovae. 
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I. INTRODUCTION 


The properties of neutrinos attract a strong attention since neutrinos alone are the el¬ 
ementary particles showing the evidence of new physics beyond the standard model. In 
particular, recent experiments of the SuperKamiokande (SK) detector [|I|] and the Sudbury 
Neutrino Observatory (SNO) [2] have shown that neutrino have non-zero mass and mixing 
angles by solar B B | and atmospheric J7|, Q neutrino data. These neutrino oscillation 
experiments enable us to constrain the squared mass difference between mass eigenstates 
and mixing angles. The solar neutrino experiments strongly support the matter-induced 
large mixing angle (LMA) solution (A/?t 2 2 = 2 — 10 x 10 -5 eV 2 , tan 2 d 12 = 0.25 — 0.5 [b]) 
and the atmospheric neutrino experiments have proved that and v T are maximally mixing 
(Am | 3 = 2 — 5 x 10 ~ 3 eV 2 , sin 2 2023 > 0.88 ||). However, there are several parameters 
which are not sufficiently constrained. One such parameter is $ 13 , for which only upper 
bound is known from reactor experiment 0 (sin 2 0 13 < 0.1). In addition, the nature of 
neutrino mass hierarchy (normal or inverted) is also matter of controversy. 

There is another astrophysical source of neutrinos, or supernova. Neutrino burst from 
the next Galactic core-collapse supernova will be captured by the present underground 
detectors, such as SK and SNO. It is widely believed that the detected events will provide 
rich information not only on the mechanism of supernovae but also on the physical properties 
of neutrinos themselves. I 11 fact, neutrino signals from SN1987A which were detected by the 
Kamiokande [llU] and IMB [O] detectors permitted many discussions about the fundamental 


properties of neutrinos |T2]. [T3], |T3)| , although the detected event numbers were quite small 
(11 for Kamiokande, 8 for IMB). What we can learn from the next Galactic supernova is 


considered in many articles (for a review, see Ref. |15]). For example, we can constrain the 
properties of neutrino oscillation such as $13 and the mass hierarchy [|H|, [17[ . 


In addition to neutrino mass, non-zero magnetic moment is another nature of neutrinos 
beyond the standard model. The most stringent upper bound of the neutrino magnetic mo¬ 
ment is obtained from cooling argument: any new physics that would enhance the neutrino 
diffusion, such that the cooling time drops below the observed duration, must be excluded. 
The recently revised upper bound is fi v < 1 - 4 x 10 ~ 12 ii b 0, where /is is the Bohr 
magneton. The above value is several orders of magnitude more stringent than the existing 
limit from reactor |L9] and solar neutrino experiments and comparable to the 
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bound from globular cluster red giant cooling arguments ||23|, which give /i v 3 x 10~ 12 /i_b. 
On the other hand, from the minimally-extended standard model calculation, we obtain 
/ij, ~ 3 x 10 - 19 //B(m 1/ /1 eV), which is many orders of magnitude smaller than the current 
experimental upper bounds. However, some particle physics models have been proposed 
to induce a larger magnetic moment of order ~ 10 _11 /i.B |23|], which is comparable to the 
current upper bounds. 

If neutrinos have the non-zero magnetic moment, it leads to precession between left- and 
right-handed neutrinos in sufficiently strong magnetic fields expected in the inner region 


of the core-collapse supernova 


]. In general, non-diagonal elements of the magnetic 


moment matrix are possible and neutrinos can be changed into different flavors and chi¬ 


ralities 127]. Furthermore, with the additional effect of the coherent forward scattering by 


matter, neutrinos can be resonantly converted into those with different chiralities |28| by a 
mechanism similar to the well known Mikheyev-Smirnov-Wolfenstein (MSW) effect |29|, |30| . 
This resonant spin-flavor (RSF) conversion induced by the neutrino magnetic moment in 
the supernova has been studied by many authors [ER |3l], |32], Fj3, [34]. |35], 


I, |38|. Almost 


all these earlier publications considered the RSF conversions between two-flavor neutrinos 
{y e v^). However, since all three-flavor neutrinos and anti-neutrinos are emitted from 
the supernova core, we should calculate conversion probabilities using three-flavor formula¬ 
tion to obtain the spectrum of the neutrinos in the detectors. Fortunately, our knowledge 
of the neutrino mixing parameters permits the three-flavor analyses, which have not been 
considered because of many unknown parameters in those days. 

In this paper, we study the RSF effect in detail, including three flavor neutrinos and 
anti-neutrinos. Particularly, a new crossing diagram, which deals with both MSW and 
RSF effects, is given. From the diagram, we find that the combination of the MSW and 
the RSF effects makes crossing schemes very interesting to investigate, and the expected 
flux would be different from that obtained in the case of pure MSW or pure RSF effect. 
In addition, we calculate numerically the conversion probabilities of neutrinos and anti¬ 
neutrinos in the supernova using three-flavor formulation. In the calculations, we adopt the 
latest neutrino mixing parameters, which are constrained by the solar and the atmospheric 
neutrino experiments, and assume the normal mass hierarchy. However, there are two 
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unknown parameters of neutrinos, /j, v and $13 1 . For / i v , we assume the value 10 ~ 12 /is 
2 , and for 6*13, two large and small values (sin 2 26*13 = 0 . 04 , 10 -6 ). Using the calculated 
conversion probabilities, we give the expected energy spectra at SK and SNO detectors, 
which show various behaviors as we change the value of p, v B. Finally, we discuss the model 
uncertainties and how to obtain the implication for the value of neutrino magnetic moment 
from the obtained energy spectra. 

This paper is organized as follows. In Section |IJ, we give the formulation used in our 
calculation, which includes all three flavor neutrinos and anti-neutrinos. In Section |TTT1 . the 
neutrino oscillation models and supernova models (original neutrino spectrum, density, Y e , 
and magnetic held profiles) considered in this paper arc illustrated and their validity is also 
discussed. In Section 0 we give a new crossing diagram, with which qualitative discussions 
concerning neutrino conversions in supernova matter become possible. The numerically 
calculated conversion probabilities and energy spectra at SK and SNO are shown in Section 
0. Finally, detailed discussions on model uncertainties and how to obtain the implication 


for neutrino magnetic moment are presented in Section VI 


II. FORMULATION 


The interaction of the magnetic moment of neutrinos and magnetic fields is described by 


((o)i? | | (l/j , ') l) 


( 1 ) 


where /i^- is the component of the neutrino magnetic moment matrix, B± the magnetic 
held transverse to the direction of propagation, (u) R and (u) L the right- and left-handed 
neutrinos, respectively, and i and j denote the havor eigenstate of neutrinos, i.e., e, //, and r. 
If neutrinos are the Dirac particles, right-handed neutrinos and left-handed anti-neutrinos 
are undetectable (sterile neutrinos), since they do not interact with matter. The spin-flavor 
conversion into these sterile neutrinos suffers strong constraints from observation of neutrinos 
from SN1987A |32|, fl4[ |39|. On the other hand, if neutrinos are the Majorana particles, 


1 There are many unknown properties in supernova physics, such as density profile, Y e , and magnetic field. 


We give in Section |I| the adopted properties in our calculation. 

2 Since the neutrino magnetic moment enters into the evolution equation only in the combination (see 
Section i>. the results would also apply to any other value of Rl/ provided that the supernova magnetic 
held strength is rescaled accordingly. 
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u R s are identical to anti-particles of u F s and interact with matter, and the constraint above 
becomes considerably weak. In this paper, we assume that neutrinos are Majorana particles. 
The diagonal magnetic moments is forbidden for the Majorana neutrinos, and therefore, only 
the conversion between different flavors is possible, e.g., (p e ) R {y^ T ) L . 

The coherent forward scattering with matter induces the effective potential for neutrinos, 
which is calculated using weak interaction theory. The effective potential due to scattering 
with electrons is given as 

V±± = ±V2Gp ^±- + 2sin 2 Ow'j n e , (2) 


where n e is electron number density, G F the Fermi coupling constant, and 9\y the Weinberg 
angle. The ± sign in front refers to v (+) and D (—) and that in the parentheses to v e (+) and 
v j) T (—). The difference between e and /i, r neutrinos comes from the existence of charged- 
current interaction. The subscript ±± of V refers to the first and the second ± sign. The 
ordinary MSW effect between v e and v^ T is caused by the potential difference, V e — V^ r = 
V ++ — = \/2G F n e . To include the RSF effect, which causes the conversion between 

neutrinos and anti-neutrinos, we should take into account the neutral-current scattering by 
nucleons: 


V = V2G f 



2sin 2 O-w^j 7i p — \/2G F -n r 


(3) 


where n p , n n are proton and neutron number density, respectively. For neutrinos we add 
+V to the potential and for anti-neutrinos —V. Therefore, the RSF conversion between u e 
and obeys the potential difference 


AR = V s - V^ T 

= (R_ + -V)- (R + _ + V) 

= V2G F -^(l-2Y e ), (4) 

m N 

where p is the density, m,v the nucleon mass, and Y e = n e /{n e + n n ) is the number of 
electrons per baryon. (When we obtain Eq. (||), we assumed charge neutrality n e = n p .) 


A. The simplest case: v e conversion 

In this subsection, we give the simplest formulation between v e and and consider the 
properties of the RSF conversion. The time evolution of the mixed state of v e and is 
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described by the Schrodinger equation 


dr 


( - 

/ 


= 

\ v n 

l 



( 5 ) 


\ "r 


where r is the radius from the center of the star, /x e/i the transition magnetic moment, and 
AH is defined by 

A "^ 2 cos 20,0 — AV. 


AH = 


2 E„ 


( 6 ) 


The resonance occurs when AH = 0, thus, the necessary condition is Y e < 0.5. (If Y e > 0.5, 
the v e <r-> P M conversion occurs conversely, at almost the same place.) Several earlier authors 
f35[, [36l pointed out that the RSF conversion occurs quite efficiently in the region above the 
iron core and below the hydrogen envelope of collapsing stars, namely, in the “isotopically 
neutral region” (O+Si, O+Ne+Mg, O+C, and He layers), in which the value of Y e is slightly 
less than 0.5 [typically, (1 — 2 Y e ) ~ 10 -4 — 10~ 3 ]. Actually, the pre-supernova model used 
in our calculation have the similar profile, as is shown in Section [111 C| . 

When the resonance is adiabatic (or when the density is slowly changing at that point 
and the magnetic held is strong enough), the significant conversion occurs. The adiabaticity 
of the RSF conversion can be derived by the same argument used in the case of the MSW 
conversion 


30, 35]. ft is given by 


rsf _ (2/i e/t -ER) 2 

/e/x 


(7) 


conversion 


conversion occurs 


\dAV/dr\ 

where the subscript “res” means that the value is evaluated at the resonance point. There¬ 
fore, if the magnetic held is sufficiently strong at the resonance point, the P, 
occurs completely. 

At the end of this subsection, we discuss the possibility that u, 
in solar magnetic held. There are several works examining the RSF conversion as a solution 
to the solar neutrino problem |4(J . In this paper, however, we postulate that the solution to 
the solar neutrino problem is given by the ordinary MSW effect (LMA solution). In fact, if 
we take the mixing parameters of LMA (A m\ 2 = 2 — 10 x 10 -5 eV 2 , tan 2 d 12 = 0.25 — 0.5) 
in calculating the RSF resonance condition AH = 0, the required density is found to be 
larger than the central density of the Sun. (In this calculation, the sign of AV is changed 
in Eq. (h;), since Y e > 0.5 in the Sun. This corresponds to u e P M conversion.) 
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B. Three-flavor formulation 


Here, we present the three-flavor (six-component) formulation of neutrino mixing by 


where 


. d 
l ~dr 



H 0 

—B i M 


B ± M 

Ho 




( v \ 


( p \ 

V = 

'A 

, v = 

A 




\ P r) 


( 8 ) 


(9) 


Hn = 


2 E v 


-U 


0 0 
0 Am^ 


0 

0 


^ (v ++ + v 


t/ f + 


Ho = 


2 E v 


-U 


y 0 0 Amyy 

/ 0 0 0 \ 
0 A m\ 2 0 


U ] + 


0 0 

0 E+_ + V 0 

0 0 V+ - + v J 

(V- + - V 0 o \ 

0 V— - V 0 



V° 

0 

Aml 3 J y 0 0 V _— V y 


1 U el u e2 u e3 ^ 


^ C12C13 S12C13 S 13 ^ 

u = 

Ufi l Ufj , 2 Ufj , 3 

= 

— S12C23 — C12S23S13 C12C23 — Sl2S23 s 13 S23 c 13 


to 

CO 


y S12S23 — C12C23S13 —C12S23 — S12C23S13 C23C13 J 


( 10 ) 


( 11 ) 


( 12 ) 


0 /^e^t fJ'er 

fJ'efi 0 i 


(13) 


y f^er fJ'fiT 0 J 

and Cij = cos 9ij, = sin 6 ir (We assume CP phase 5 = 0 in Eq. ( |T2| ) for simplicity.) 

The resonant flavor conversion occurs when the two diagonal elements in matrix in Eq. (^J) 
have the same value. There are five such resonance points, which are for v e (MSW- 

L), v e v r (MSW-H), v e <-► Up (RSF-L), v e <-► v T (RSF-H), and v T conversions. 

Suffixes “-L” and “-H” attached to “MSW” and “RSF” indicate whether the density at the 
resonance points are lower or higher. Hereafter, we neglect <->■ u T conversion, since it is 
always nonadiabatic as well as including this makes discussion further complicated. 
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III. MODELS 


A. Neutrino parameters 


We adopt the neutrino mixing parameters introduced to explain the observations of the 
solar and the atmospheric neutrino experiments i i a a a a, and assume the normal 
mass hierarchy. From the atmospheric neutrino experiments, we use the values Q 

Am 2 l3 = 2.8 x 10" 3 eV 2 , sin 2 20 23 = 1.0, (14) 

and from the solar neutrino experiments ||, 

Am 2 2 = 5.0 x HP 5 eV 2 , tan 2 d l2 = 0.34, (15) 


in our calculations. For 0 13 , which is not sufficiently constrained, we adopt two large and 
small values, 

sin 2 2013 = 0.04, 1.0 x 10" 6 , (16) 


where each value causes the adiabatic and nonadiabatic MSW-H conversions, respectively. 

As for neutrino magnetic moment, there are three transition moments for Majorana 
neutrino (Eq. ©)■ We assume that all of them are near the present upper bound (i.e., 
/j,ij = 10 ~ 12 hb)- However, as we will show in Section [111 D| , /r eT alone governs results. 


B. Original spectrum of neutrinos 


We use a realistic model of a collapse-driven supernova by the Lawrence Livermore group 
El to calculate the neutrino energy spectrum at production. We show in Fig. |T] the 
time-integrated energy spectra of neutrinos (see Ref. for detail), where the distance 
to the supernova is set to 10 kpc. (For discussions from here, we adopt this assumption, 
D — 10 kpc.) Both the flux integrated over entire time during neutrino burst and that 
integrated until 0.5 seconds after core bounce are shown in the same figure. Although the 
neutrino burst persists ~ 10 seconds, about half of neutrinos are emitted by 0.5 seconds, 
according to numerical simulation by Lawrence Livermore group [Q. The reason for using 
time integrated flux until 0.5 seconds is given in the next subsection. 
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The mean energies are different between flavors: 

{E Ue ) ~ 11 MeV, (17) 

(E Pe ) ~ 16 MeV, (18) 

(E Ux ) ~ 22 MeV, (19) 

for fully time-integrated spectrum, while for that integrated until 0.5 seconds, 

(E Ve ) ~ 11 MeV, (20) 

(Ep B ) ~ 14 MeV, (21) 

(TU ^ 18 MeV. (22) 


This hierarchy of mean energies is explained as follows. Since z/^’s interact with matter only 
through the neutral-current reactions in supernova, they are weakly coupled with matter 
compared to yds and z? e ’s. Thus the neutrino sphere of z/^’s is deeper in the core than that 
of z/ e ’s and z/ e ’s, which leads to higher temperatures for u^s. The difference between z/ e ’s and 
P e ’s comes from that the core is neutron rich and z/ e ’s couples with matter more strongly, 
through u f /ri —* e~p reactions. 


C. Supernova progenitor model 


We use the precollapse model of massive stars of Woosley and Weaver [|7| . The progenitor 
mass was set to be 15M 0 , and the metalicity the same as that of the Sun. The upper panel of 
Fig. |2] shows pY e and p{l — 2Y e ) profiles of this progenitor model. In the same figure, we also 
show A 12 = Am( 2 cos 20i 2 /2^2GfE u and Ai 3 = Amf 3 cos2#i 3 /2\/2G(Fi7 i ,. At intersections 
between Ai 2 , A i3 and p( 1 — 2Y e ), pY e , the RSF and MSW conversions occur (see Eq. (|6])). 

However, in fact, since the density profile changes drastically during neutrino burst (~ 10 
sec) owing to shock propagation, we cannot use the static progenitor model when calculating 
the conversion probability in supernova ^4j. We show in Fig. § the time-dependent p(l—2Y e ) 
and pY e profile calculated by Lawrence Livermore group |45|, compared with the static model 
adopted in our calculation [43|. From the figure, the shock wave changes the p(l—2V e ) profile 
in the resonance region, which is responsible for RSF conversions, at > 0.5 seconds after 
bounce, while that at < 0.5 seconds is almost the same as the static progenitor model. Thus, 
from here, we confine our discussion until 0.5 seconds after core bounce, since in that case, 
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using the static progenitor model is considered to be a good approximation. Calculating for 
full time-scale during neutrino burst with the time-dependent density and Y e profiles makes 
discussion very complicated, and it is out of scope of this study, while it will be our future 
work. 


D. Supernova magnetic field 


We assume that the global structure of the magnetic field is a dipole moment and the 
field strength is normalized at the surface of the iron core with the value of 10 8 — 10 10 
Gauss. The reason for this normalization is as follows. As pointed out by Totani and Sato 
36| , the global magnetic field normalized at the surface of the newly born neutron star, 


such as used in Ref. [[I|], is not a good assumption. It is because that the magnetic fields 
of a nascent neutron star are hard to affect the far outer region, such as the isotopically 
neutral region, within the short time scale of neutrino burst. The magnetic fields, therefore, 
should be normalized by the fields which are static and existent before the core collapse. In 
addition, as discussed in the previous subsection, since the shock wave does not affect the 
resonance region at 0.5 seconds after bounce, it is also expected that the magnetic field 
strength at the resonance points is not seriously changed at that time. The strength of such 
magnetic fields above the surface of the iron core may be inferred from the observation of 
the surface of white dwarfs, since both are sustained against the gravitational collapse by 
the degenerate pressure of electrons. The observations of the magnetic fields in white dwarfs 
show that the strength spreads in a wide range of 10' — 10 9 Gauss |46|]. Taking account of 
the possibility of the decay of magnetic fields in white dwarfs, it is not unnatural to consider 
the magnetic fields up to 10 10 Gauss at the surface of the iron core. Then, in Eq. (Q, 
Bj_ = B 0 (r 0 /r) 3 sin 0, where B 0 is the strength of the magnetic field at the equator on the 
iron core surface, r 0 the radius of the iron core, and 0 the angle between the pole of the 
magnetic dipole and the direction of neutrino propagation. Hereafter, we assume sin© = 1. 

From the magnetic field, density, and Y e profiles of supernova we can calculate the adia- 
baticity parameter of each resonance points with Eq. (|7|). The lower panel of Fig. |j shows 
adiabaticity parameters of RSF conversions as a function of radius. Those for both Bo = 10 8 
and Bq = 10 10 Gauss, when /i u = lCT 12 /^ is assumed, are shown. From the figure, we see 
that the RSF-H conversion can be adiabatic when the B 0 ~ 10 10 G, while the RSF-L conver- 
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sion is always nonadiabatic under any conditions considered here. Therefore, the values of 
transition magnetic moment yU e/i , /i MT do not matter, while that of fi eT alone governs results. 


IV. LEVEL CROSSING SCHEMES AND QUALITATIVE ILLUSTRATIONS 


From the formulation which includes three-flavor neutrinos and anti-neutrinos (Eq. (||)), 
we give a new level crossing diagram, which is an expanded version of the ordinary diagram 
including MSW resonance alone (see Fig. 5 in Ref. |L 6 |). We show it in Fig. f|, in which 
the adiabatic conversion means that the neutrinos trace the solid curve at each resonance 
point (i.e., the mass eigenstate does not flip), while the nonadiabatic conversion the dotted 
line. The figure clearly includes not only ordinary MSW resonances but also the RSF 
effects, and it is expected that the combination effect of MSW and RSF makes this scheme 
very interesting to investigate. From the earlier publications on the MSW effect, we know 
that the MSW-L is adiabatic for LMA solution, and MSW-H is adiabatic for large $13 and 

T7fl. The adiabaticity of the RSF-L and RSF-H depends on 


nonadiabatic for small $13 


the magnetic field strength and the density profile at the resonance point (see Eq. (0)), 
which are given in Section |TTI] and Fig. 0. 

Here, we qualitatively consider the flavor transition of neutrinos under some simple con¬ 
ditions. At the supernova core, electron neutrinos and anti-neutrinos are produced as the 
mass eigenstate in matter, owing to large matter potential. The other flavor (/x, r) neutrinos 
and anti-neutrinos at production are not mass eigenstates, however, we can make them by 
linear combination of these flavor eigenstates. (These states corresponds to z/, z/, P f , v' T in 
Fig. |.) Then, neutrinos and anti-neutrinos of all flavors propagate in supernova matter 
along the curves in Fig. [4j. being affected by four resonances, and appear from the stellar 
surface as pure vacuum mass eigenstates or their mixing states. 

When both the RSF-L and RSF-H are nonadiabatic, the situation is the same as the pure 
MSW case, i.e., the conversions occurs as 


*4 *4 


"'a -»■ *4 K 


V2 


(adiabatic MSW-H), 




Rz *4 


*4 
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Ve 


V 2 V e 
V 1 ^ 
^3 K 


V\ 

V2 

"3 


(nonadiabatic MSW-H). 


From these considerations, we can estimate the spectral deformation of v e and v e , which are 
the easiest to detect at SK and SNO. When the MSW-ff is adiabatic, we obtain 

F e = \U el \ 2 F 1 + \U e 2 \ 2 F 2 +\U a \ 2 F 3 
= + \U a \ 2 Ft + \U ei \ 2 F^ 

= (1 - \Ua\*)F? + \Ua\F?, (23) 

F, = |tr«i| 2 it+(l-|tf.i| 2 )fj, (24) 

where F, is the flux of the i-type neutrinos at detection, F® the flux at production which is 
shown in Fig. |I|. and subscript x — n,r, £l, f, since the flux of these flavors are considered 
to be the same. With the same manner, for nonadiabatic MSW-ff, we obtain 


F e = (l.»\U e2 \ 2 )F°+\U e2 \ 2 F°, 
F- e = \U e i\ 2 Fg + (1 — \U e i\ 2 )F® 


(25) 

(26) 


When the RSF-ff is adiabatic and RSF-L nonadiabatic, which is actually the case for 
strong magnetic field (see Fig. ||), the same arguments lead to for adiabatic MSW-ff: 


v? 


z/ 3 v e —> v 2 
vi v' tl v 2 
V\ K —> V 3 


F e = \U el \ 2 F° + \U e2 \ 2 F" + \U ( 


2 p 0 


F- = F 

r e r x ? 


and for nonadiabatic MSW-ff: 




e3 


| 2 F ° 


V2 V e -> V 3 
-! 


Vi w -> ^ 2 
z?' -> i/ 3 


(27) 

(28) 


F e = \Uei\ 2 F° x + |f/ e2 | 2 F e u + iFesI 2 ^, 


2 77 O 


1 2 77O 


Fe 


= 


(29) 

(30) 
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The adiabatic RSF-H conversion, thus, deforms the energy spectra, and the difference from 
the case of the pure MSW effect is clear. 

V. RESULTS OF NUMERICAL CALCULATIONS 

A. Conversion probability 

We calculated Eq. (0) numerically with adopted models given in Section ES and obtained 
the conversion probabilities for each flavor. Figure 0 shows those for P e (upper panel (a)) 
and Up (lower panel (b)), under the condition that B 0 = 4 x 10 9 G, E u = 20 MeV, and 
sin 2 2013 = 0.04. Since z/ T ’s (h T ’s) are maximally mixing with v^s (P M ’s), we did not show 
P( v e —> z/ r ,P T ) and P(i/ M —> u T ,u T ), which are almost the same as that for We 

can see three clear conversions, one RSF and two MSWs. The RSF-ff conversion occurs 
at r ~ O.OO7P 0 between P e and v T . (The apparent P e conversion occurs, since /y t 's 

are maximally mixing with z/ r ’s.) After that, neutrinos propagate in matter inducing the 
ordinary MSW effects, which are studied in detail by earlier authors |TB], [T7| (MSW-H at 
~ O.O5P 0 , MSW-L at ~ (0.2 — 0.5)72©). The strange behavior of P(r / (U —> P e ) at r ~ O.O4P 0 
is considered due to the nonadiabatic RSF-L conversion, while the probability is not changed 
as a whole. 

Figure | 6 ] is the same one as Fig. 0, but it is for sin 2 29i 3 = 10 -6 instead of 0.04. In this 
case, we do not see the MSW-H conversion, since it is nonadiabatic. 

We show in Fig. |7j. the dependence of the conversion probability P(z/ ;i —> u e ) on the 
energy and the magnetic field strength at the iron core, assuming sin 2 2^13 = 0.04. From 
Fig- 0 (a), we see that the resonance point goes to outer region and the conversion becomes 
more and more nonadiabatic, as the neutrino energy becomes larger. In the same manner, 
Fig. 0 (b) indicates that the conversion becomes more adiabatic for larger magnetic field 
strength. These behaviors can be understood qualitatively using Fig. j2| and Eq. (0)- 

At the end of this subsection, we compare the calculated conversion probabilities with 
ones which are evaluated with time-dependent density profile. Figure 0 shows P(P e —> P e ) 
and P(P e — 1 ► v^r) as a function of neutrino energy. Since how the shock propagation changes 
the structure of magnetic field at resonance points is not known well, we assume the same 
static dipole profile as for magnetic held, and the figure is for B 0 = 10 10 G. As discussed in 
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Section [Til D| . this assumption is considered to be a good one for < 0.5 seconds after bounce, 
however, probably incorrect for > 0.5 seconds. Anyway, from the figure, we can confirm that 
the probabilities are in good agreement between the model 0.5 sec and the static progenitor 
model. 


B. Energy spectrum at the SuperKamiokande detector 


SK is a water Cherenkov detector with 32,000 tons of pure water based at Kamioka in 
Japan. The relevant interactions of neutrinos with water are as follows: 


& e +P ~ 

-> e + 

+ n, 

(C.C.) 

(31) 

v e + e ~ - 

+ "e 

+ e“, 

(C.C. and N.C.) 

(32) 

v e + e ~ - 

Ve 

+ e j 

(C.C. and N.C.) 

(33) 

u x + e ~ - 

v x 

+ e", 

(N.C.) 

(34) 

v e + o - 

e~ 

+ F, 

(C.C.) 

(35) 

v e + o - 

+ e + 

+ N, 

(C.C.) 

(36) 


where C.C. and N.C. stand for charged- and neutral-current interactions, respectively. 

SK is not in operation now because of the unfortunate accident, which occurred on 12th 
November 2001. However, it is expected to be repaired using remaining photomultiplier 
tubes (PMTs), and will restart the detection by the end of this year (2002) with lower 
performance. The effect of the accident on its performance is expected not to be serious 
for supernova neutrinos, because the fiducial volume will not change, and the threshold 
energy change (from 5 MeV to about 7-8 MeV) influences the event number very little. 
Although the energy resolution will become about \/2 times worse, it does not matter for our 
consideration. In the calculations, we used the energy threshold and the energy resolution 
after the accident. 

For the reaction cross sections, we refer to Ref. [[Tj. The v e p is the largest reaction 
and the energy of recoil positrons is related to that of injected neutrinos by simple formula 
E e = E Pe — ( M n — M p ) = Ep e — 1.29 MeV, where M n , M p are neutron and proton mass, 
respectively. The energy spectrum of electrons (positrons), therefore, reflects that of u e , 
which is easily estimated in some simple cases using e.g., Eqs. ([24]), ([261), ([281), and fl30[). 
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Figure [| shows the energy spectra of electrons (positrons) for first 0.5 seconds evaluated 
with sin 2 2013 = 0.04 and several values of Bq, which are calculated with conversion proba¬ 
bilities, original neutrino spectra, reaction cross sections, and detector performance. When 
we increase the value of B 0 , the energy spectrum shifts to higher energy region and the shift 
saturates at B 0 ~ 10 10 G. This is because at B 0 > 10 9 G the RSF-H conversion starts 
to become adiabatic and it achieves completely adiabatic transition at Bo rv./ 10 10 G. For 
sin 2 2013 = 10” 6 , we obtained almost the same figure as Fig. |9|. and we do not show it here. 


C. Energy spectrum at the Sudbury Neutrino Observatory 


SNO is a Cherenkov detector filled with 1,000 tons of heavy water. The interactions of 
neutrinos with heavy water are 


v e + d - 

■* e +p + p, (C.C.) 

(37) 

v e + d - 

+ e + + n + n, (C.C.) 

(38) 

v x + d - 

■+ v x +p + n, (N.C.) 

(39) 


and reactions represented by Eqs. (^) to (|36|). Furthermore, the produced neutrons are 


detected through the interaction with surrounding nuclei as 

7i + d, —> 3 H + 7 6i3 , (efficiency: 24%) 
n + 35 C1 —■> 36 C1 + 7 8 . 6 , (efficiency: 83%) 


(40) 

(41) 


where 76.3 and 73.6 represent the 7 -rays whose energies are 6.3 MeV and 8.6 MeV, respec¬ 
tively, and NaCl is added to efficiently capture neutrons. Owing to these neutron capture 
reactions, the neutrino detection through the neutral-current reaction Eq. (|39|) is possible. 
Furthermore from the charged-current reaction of z? e ’s (Eq. (|38D), we receive more than two 
signals, i.e., one from the preceding recoil positron, and the others from the delayed 7 -rays 
by neutron captures. Thus, using this criterion, we can discern the z7 e d events from the other 
detection channels. 

For the cross sections of reactions with deuterons, we referred to Ref. [47]. Figure IT 
shows the energy spectra for first 0.5 seconds evaluated with sin 2 2013 = 0.04 and several 
values of B 0 3 . The upper panel (a) shows the total energy spectra and the lower panel 


Since the neutral-current reaction Eq. (]39|) does not contribute to energy spectrum, we neglect it from 
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(b) the spectra of (total— v e d) events which can be obtained using the delayed coincidence 
signals, where the line types are the same as those used in Fig. As Bq inceases, the energy 
spectrum of the total events becomes harder, while that of the (total —D e d) events softer. 
These properties can be explained as follows. Since almost all the events come from v e d and 
v e d reactions, the (total— u e d) events are roughly coincide with v e d events. Thus, from the 
qualitative discussion given in Section B we can obtain the properties in Fig. H (b) (see 
Eqs. (^3|) and (|7|)). Since the energy spectrum due to v e d becomes harder for larger B 0 
(see Eqs. (^4|) and (p8|)), the total spectrum also becomes harder but with the suppressed 


degree of the shift, because that of the u e d events relaxes. 

Figure [Q] is the same one as Fig. 0 but for the different value of 013 (sill 2 26» 13 = 10 6 ). 
For larger B 0 , the spectrum of the total events becomes harder, while that of (total— D e d) 


events changes little (see Eqs. (|25|), (p6|), fl29|), and (p0|)). 


VI. DISCUSSION 


In this study, we used the original neutrino spectrum calculated by Lawrence Livermore 
group |42|]. Unfortunately, their study as well as the other published full numerical supernova 
collapse simulations have not yet included the nucleon bremsstrahlung process or nucleon 
recoils, even though it is no longer controversial that these effects are important. Recent 
studies (e.g. J48|) including all these processes have shown that average v x energy exceeds 
the average v e energy by only a small amount, 10% being a typical number. If this is the case, 
the spectral deformation, particularly at SK detector, will be little and it is difficult to obtain 
some information on the flavor conversion mechanism in supernova. However, it is premature 
to conclude that their results are correct, since it is based on the neutrino transport study 
on the background of an assumed neutron star atmosphere, and this approach lacks hydro- 
dynamical self-consistency. Further it is also because the mean energies and their ratios 
change significantly between the supernova bounce, accretion phase, and the later neutron 
star cooling phase. 

Thus, at present, even if the Galactic supernova neutrino burst were detected by SK and 
SNO, it would be very difficult to obtain useful information on whether the RSF effect actu- 


here on. 
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ally occurs in supernova and to obtain implication for neutrino magnetic moment, owing to 
large uncertainties concerning original neutrino spectrum. In the future, however, it might 
be plausible that these uncertainties will be removed by development of numerical simula¬ 
tions, and the signal of supernova neutrino burst would provide valuable implications for 
the flavor conversion mechanisms. In discussions from here, we present the method to probe 
conversion mechanisms from the detected neutrino signals, assuming that the uncertainties 
concerning the supernova models are considerably reduced by the future development of 
numerical approaches. 


A. Implication for neutrino magnetic moment 


As we have shown in Section [V], the deformation of spectral shape occurs when B 0 > 
10 9 G. Therefore the ratio of high energy events to low energy ones is a good measure 
to investigate the magnetic field strength in supernova or neutrino magnetic moment. We 
adopted as that two quantities below: 

Number of events for E P > 25 MeV 


Rsk — 


R. 


Number of events for E e < 20 MeV 
Number of events for E e > 25 MeV 


SNO 


(42) 

Number of events for E e < 15 MeV ^ ^ 

Figure |12] shows Ask as a function of B 0 sin 0 for sin 2 20 13 = 0.04, where error bars include 
only statistical errors. (Here, we re-introduce the angular factor sin 0 since we also discuss 
the dependence on magnetic field orientation below.) We obtained almost the same figure 
for sin 2 26 i 3 = 10“ 6 . Therefore, SK data will provide statistically sufficient information 
on the magnetic held strength in the supernova and neutrino magnetic moment, and these 
consequences are independent of ffi 3 . For example, the value of Ask which is found to be 
larger than ~ 1 indicates that the RSF conversion occurred in the supernova, and we can 
constrain the value of fi u B 0 sin 0 assuming the supernova and presupernova models. 

The same figure for SNO is Fig. 13, where the upper panel (a) is for sin 2 20 13 = 0.04 
and the lower panel (b) for sin 2 2ffi 3 = 10~ 6 . If we obtain Asno (total) which is larger than 
Rsno( total — z/ e d), it suggests that the supernova magnetic held is quite strong, B 0 sin© > 
5 x 10 9 G, and neutrinos have non-zero magnetic moment. However, since the statistical 
errors of Rsno are rather large because of smallness of data, we can use the SNO data to 
conhrm the SK result. 
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In addition, if we know the value of the angular factor sin 0 from the optical observation 
such as the Hubble Space Telescope (HST), we can constrain the magnetic held strength 
Bq (not the combination B 0 sin 0) from the neutrino observation. Actually, in the case of 
SN1987A, which occurred in the Large Magellanic Cloud (LMC), 0 is constrained with the 
procedure as follows. The HST images of SN1987A show a bright elliptical ring surrounding 
the optical remnant. This inner ring has been interpreted as a ring of high-density material 
in the equatorial plane of the progenitor, caused by the impact of the high velocity stellar 
wind of the blue supergiant progenitor with its earlier, slower red giant wind, and ionized 
by the UV flash from the supernova explosion j49[. The inner ring is known to be inclined 
to the line of sight by ~ 43°. Thus, assuming that the magnetic dipole moment of the star 
orient to the rotation axis, we can constrain sin© to be ~ 1/a/2 from the HST observation. 
However, if the next supernova occurred in the Galactic center, it would be quite difficult to 
estimate 0 since the Galactic center is optically thick from the Earth. On the other hand, 
if it occurred in the optically thin environment such as LMC, we would be able to constrain 
0 from optical observations and Bq from the neutrino signals. Although LMC is located 
at ~ 50 kpc away from the Earth, SK data would permit statistically sufficient discussions 
on the supernova magnetic field. For example, in that case, the la statistical errors of i?gK 
are as large as 5a errors in the case of the Galactic supernova (D = 10 kpc), which are also 
shown in Fig. [T2|, and it will be sufficient to distinguish whether 7? o sin0 > 5 x 10 9 G or 
Bq sin 0 < 10 9 G. 


B. Constraint from supernova relic neutrino observation 


In addition to next Galactic supernova, the observation of supernova relic neutrino (SRN) 
background would be available to obtain information on flavor conversion mechanism in 
supernova. The flux of SRN is calculated by 

dF v r 


r max d ( \ dNu(E' ) dt 

dE= C i + 


(44) 


where E' v = (1 + z)E v , Rs^(z) is supernova rate per comoving volume at redshift z, dN v /dE v 
energy spectrum of emitted neutrinos, and 2 max the redshift when the gravitational collapses 
began (see Ref. |)0) for detailed discussion). 

Most recently, SK collaboration set an upper bound of 1.2 P e cnr 2 s -1 for the SRN flux 
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in the energy region E u > 19.3 MeV |5T|. That limit is the same order as some typical 
theoretical predictions such as Ref. |50|. For example, Ando et al. |5(J predicted that the 
total SRN flux integrated over entire energy was 11 cnr 2 s _1 , while the corresponding SK 
limit calculated with the predicted spectral shape is 31 cm -2 s -1 . Since the theoretical cal¬ 
culations contain many ambiguities such as the supernova rate in the universe and neutrino 
spectrum from each supernova, this severe observational SRN limit can provide a number 
of valuable information on the various fields of astrophysics and cosmology. Further, in the 
near future, it is expected that the upper limit will be much lower (about factor 3) when the 
current SK data of 1,496 days are reanalyzed using some technique to reduce background 


against detection |52|]. In that case, the further severer constraint not only on the star 
formation history but also on the nature of neutrinos might be obtained. 

Since the RSF mechanism further enhances the average energy of u e compared to the 
ordinary MSW effect, the current and the future expected upper limit on the SRN flux 
might provide very useful information on that mechanism. However, in order to give precise 
SRN prediction including the RSF mechanism, we must calculate conversion probability in 
supernova for full time-scale during the neutrino burst. Therefore, the prediction is difficult 
at present stage, and it is also one of our future works. 


VII. CONCLUSION 

We investigated three-flavor neutrino oscillation with non-zero neutrino magnetic moment 
in supernova magnetic fields. We gave a new crossing diagram which includes not only 
ordinary MSW resonance but also magnetically-induced RSF effect. From the diagram, it 
was found that there occurred four resonances, two were induced by the RSF effect and 
two by the pure MSW. In addition, the qualitative behavior of neutrino conversions and 
expected neutrino spectrum at the detector were found to be well illustrated. 

We also numerically calculated conversion probability in the supernova to obtain the 
expected energy spectrum at the detector. In that calculation, realistic neutrino mixing 
parameters inferred by recent neutrino oscillation experiments and realistic density and Y e 
profiles of presupernova stars were adopted. As for magnetic held in the supernova, we 
assumed the dipole structure, and that the held strength was normalized at the surface 
of the iron core, which was inferred from the observations of white dwarfs. Although the 
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shock wave propagation changes density profile in the resonance region drastically and is 
also expected to change magnetic hied profile at the same region, we showed that if we 
confine onr discussion until 0.5 seconds after core bounce, using static progenitor model and 
static magnetic held structure is expected to be a good approximation. 

As a result of numerical calculations, it was found that if the magnetic held at the surface 
of the iron core is sufficiently strong (B 0 ~ 10 10 G, where we assumed /i u = 10^ 12 /ib), the 
RSF-H conversion becomes adiabatic and complete V e u T conversion occurs. Hence, non¬ 
zero neutrino magnetic moment with the strong magnetic held enhances the energies of u e 
and hardens the energy spectrum at SK and SNO. Thus, in the future, if the uncertainty 
concerning original neutrino spectrum is removed by development of numerical simulations 
for supernova explosions, we might be able to obtain information on neutrino magnetic 
moment as well as the magnetic held strength in supernova using the energy spectrum at 
detectors. In addition, it could be also plausible using the current and the future supernova 
relic neutrino observations. 
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FIG. 1: Original time-integrated energy spectra of v e (solid curve), u e (dashed curve), and v x 
(dotted curve), where u x represents ;y /t . r and v llT . For each line type, the larger one represents 
fully time-integrated spectrum, while smaller one spectrum integrated until 0.5 seconds after core 
bounce. 


23 








0.001 


0.01 


0.1 

Radius 



10 


FIG. 2: Supernova profiles used in our calculations |43|]. Upper panel: The solid curve repre¬ 
sents the density and Y e combination which is responsible for the RSF conversions (p( 1 — 2Y' e )), 
and the dashed curve for the MSW conversions ( pY e ). Two horizontal bands are for A 12 = 
Amf 2 cos 2612 / 2 \/ 2 GfE u and A 13 = Amf 3 cos 28 \ 3 / 2 ^ 2 GfE u . (The band width comes from 
the energy range 5-70 MeV.) At intersections between the solid, dashed curves and the horizontal 
bands, the RSF and MSW conversions occur (see Eq. (jtj)). Lower panel: The adiabaticity param¬ 
eter of RSF conversions (Eq. ([?])) as a function of radius, where the indicated values are magnetic 
field strength at the iron core ( Bq ) when p u = 10 ~ 12 pb is assumed. 
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FIG. 3: Time-dependent density and Y e profiles of supernova calculated by Lawrence Livermore 
group |45|]. Upper panel: The p{ 1 — 2 Y e ) profile of supernova, which is responsible for RSF con¬ 
versions. The profile at 0.5, 0.7, 2, and 10 seconds after bounce, are shown. The profile calculated 
with static progenitor model [^3| is also shown, for comparison. The resonance points, which are 
represented by the intersections between the p{ 1 — 2 Y e ) curves and horizontal bands (A 13 , A 12 ), 
are affected by the shock propagation at > 0.5 seconds. Lower panel: The same figure, but for pY e 
profile, which is responsible for MSW conversions. 
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FIG. 4 : The schematic illustration of level crossings, where ^1,2,3 and Pi 2 3 represent the mass 
eigenstates of neutrinos and anti-neutrinos in matter, respectively, and v' IJ T and v' the mass 
eigenstates at production, which are superpositions of and u T or and P T . There are four 
resonance points, MSW-L, MSW-H, RSF-L, and RSF-H. The adiabatic conversion means that the 
neutrinos trace the solid curve at each resonance point (i.e., the mass eigenstate does not flip), 
while the nonadiabatic conversion the dotted line. 
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FIG. 5: (a) Conversion probabilities of v e under the condition that Bq = 4 x 10 9 G, E v = 20 MeV, 
and sin 2 2013 = 0.04. The solid curve represents P(v e —> v e ), the dotted curve P(v e —> the 
short-dashed curve P(v e —> 9 e ), and the long-dashed curve P(u e —> 9^). (b) The same figure as 
(a), but for those of The solid curve represents P(z/ M —» u e ), the dotted curve Piy^ —> v^), the 
short-dashed curve P(r / At —> 9 e ), and the long-dashed curve P(r' /i —> 9^). 
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FIG. 6 : The same figure as Fig. ||, but for sin 2 26*13 


10 -6 . 



















0.01 0.1 1 

Radius [R Q 


FIG. 7: Conversion probability P —► v e ) for sin 2 2^13 = 0.04. (a) The energy dependence 
for Bq = 2 x 10 9 G. (b) The dependence on the magnetic field strength at the iron core, for 


E u = 25 MeV. 
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FIG. 8: The conversion probability P(P e —► P e ) and P(P e —> u h,t) as a function of neutrino energy, 
which are evaluated with time-dependent density and Y e profiles at 0.5 (solid line), 0.6 (short- 
dashed line), 0.7 (dotted line), and 2 (long-dashed line) seconds after bounce. Those obtained with 
static progenitor model is also shown as dot-dashed line. The magnetic field profile is assumed to 
be the same one as the static dipole model (Bq = 10 10 G). 
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FIG. 9: Energy spectra of electrons (positrons) at SK for first 0.5 seconds, evaluated with 
sin 2 2013 = 0-04 and several values of Bq. The energy spectra for sin 2 2013 = 10~ 6 are almost 
the same. The mean energy increases as Bq becomes larger, and saturates around Bq = 10 10 G 
(completely adiabatic RSF-H conversion). 


31 







0 


10 


20 


30 


E [MeV; 


40 


50 


60 


FIG. 10: Energy spectra at SNO for first 0.5 seconds evaluated with sin 2 2f?i3 = 0.04 and several 
values of B$. The line types are the same as those used in Fig. |9[ (a) Total numbers of events 
at SNO. As £>o increases, the spectrum becomes harder, (b) Numbers of events at SNO, where 
the v e d events (Eq. (|38|)) are subtracted from the total events. The spectrum moves conversely 
compared to the total events (a), as the values of Bq. 
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FIG. 11: The same figure as Fig. 0 but for sin 2 2$i3 = 10 6 . For larger Bq, the spectrum of the 
total events (a) becomes harder, while that of (total— v e d) events (b) changes little. 
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FIG. 12: Ratio of high energy events to low energy events at SK, Rsk, as a function of Rosin© 
for sin 2 2#i3 = 0.04, where error bars include only statistical errors. The solid error bars are at 1 <t 
level and the dashed ones at 5 a level. The thick and thin bands are the cubic spline interpolations 
of la and 5 a error bars, respectively. 
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FIG. 13: Ratio of high energy events to low energy events at SNO, -Rsno> where the upper panel 
(a) is for sin 2 26U = 0.04 and the lower panel (b) for sin 2 2$i3 = 10 -6 . The error bars include only 
statistical errors, and are at la level. The thick and thin bands are the cubic spline interpolations 
of these error bars. 



35 








































